Brownian Motion on a Sphere: Distribution of Solid Angles 
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Abstract 



We study the diffusion of Brownian particles on the surface of a sphere and compute the 
distribution of solid angles enclosed by the diffusing particles. This function describes 
the distribution of geometric phases in two state quantum systems (or polarised light) 
undergoing random evolution. Our results are also relevant to recent experiments which 
observe the Brownian motion of molecules on curved surfaces like micelles and biological 
membranes. Our theoretical analysis agrees well with the results of computer experiments. 
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Let a diffusing particle start from the north pole of a sphere at time r = 0. We 
join the final position of the particle at time (3 to its initial one (the north pole) by the 
shorter geodesic. This rule is well defined, unless the final position is exactly at the south 
pole, a zero probability event. The path followed by the diffusing particle (closed by the 
geodesic rule) encloses a solid angle 0. The question we address is: At time /3, what is 
the distribution P@(Q) of solid angles? 

An experimental motivation for this question comes from recent time-resolved fluores- 
cence studies |]l|, |2| ||, U on the Brownian motion of rod- like molecules on curved surfaces 
such as micelles and lipid vesicles. The experimentally measured fluorescence anisotropy 
is affected by the curvature of the surface due to geometric phase effects fl], [|, [|. The 
problem of diffusion on curved surfaces reappears in Nuclear Magnetic Resonance (NMR) 
and Electron Spin Resonance (ESR) studies on micelles and lipid vesicles. The curvature 
of the surface affects the spin relaxation times jj], ||, |||. In biological systems, the trans- 
lational diffusion of solutes bound to various curved surfaces influences metabolic rates 
and transmission rates of chemical signals Jl], ||]. Recently, fluorescence anisotropy decay 
of molecules on curved surfaces has been studied theoretically and using Monte Carlo 
simulations (lj. This study shows that the geometric phase crucially affects the anisotropy 
decay if the molecules are tilted relative to the surface normal. We therefore need to the- 
oretically evaluate the distribution of geometric phases or equivalently, the distribution of 
solid angles in order to understand diffusion processes on curved surfaces. 

In this paper we present an answer to the question posed above and compare our 
theoretical results with Monte Carlo simulations. The theoretical analysis presented here 



is based on earlier work 11 , [T^, 13 , |Tj] . In pfl , [Lfl the distribution of solid angles for 



closed Brownian paths on the sphere is computed. However, closed Brownian paths are 
a set of measure zero among Brownian paths and a comparison of the results of fi"l| , |l~4|| 
with real or computer experiments is hampered by poor statistics. We adapt the methods 



of [11, [L4| to allow for open Brownian paths and compute the distribution of solid angles 
for such paths, closing the path by the geodesic rule || in order to make the solid angle 
a well-defined quantity. The qualitative idea behind |l0| , 11, 12, 13, [l4j] and the present 
paper is to use a magnetic field as a "counter" , to measure the solid angle enclosed in a 
Brownian Motion. 

We first illustrate our method by solving a similar problem on the plane. Let a diffusing 
particle start from the origin r = of the plane at time r = 0. Let us suppose that the 
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particle arrives at ft at time (3. Join fj to the origin by a straight line. The open Brownian 
path r(r) closed by a straight line encloses an area A. What is the probability distribution 
P(A) of areas? By "area" we mean the algebraic area, including sign. Area enclosed to 
the left of the diffusing particle counts as positive and area to the right as negative. This 
problem has been posed and solved earlier. We illustrate our method by solving this 
problem before moving to the problem of main interest in this paper. 

Let {r(r) = {x(r), y(r)}, < r < /3 , r(0) = 0} be a realization of a Brownian path 
on the plane which starts at the origin. As is well known, Brownian paths are distributed 
according to the Wiener Measure if J-[f(T)] is any functional on paths, the expectation 
value of J- is given by 

{df f S W(T)]expl-fg{i%.%dT}} 

In Eq. (|l]) the symbol 2?[r(r)] denotes a functional integral |l7j over all paths which 
start at the origin and end at ff. Finally, the endpoint ff is also integrated over. (We 
set the diffusion constant equal to half.) (3 is the time for which the diffusion has occured 
[jl8|| . Let _4[r(r)] be the algebraic area enclosed by the path r(r). Clearly, the normalized 
probability distribution of areas P(A) is given by 

P(A):=<6(A[?(t)}-A)> w . (2) 

The expectation value 4> of any function 4>(A) of the area is given by J P(A)<p(A)dA. As is 
usual in probability theory, we focus on the generating function P(B) of the distribution 
P(A) : 



P{B) := e iBA = J P(A)e l dA, (3) 

which is simply the Fourier transform of P(A). For reasons that will soon be clear, we use 
the symbol B for the Fourier transform variable. The distribution P(A) can be recovered 
from its generating function by an inverse Fourier transform. From Eqs.(j2|) and (||) above 
we find 

P{B) =< e iBA > w (4) 

Let us introduce a fictitious vector potential A.dr = (B /2)(xdy — ydx). Clearly 
§ A.dr = BA, where the integral is over the closed circuit consisting of the Brownian 
path closed by the geodesic. A has been chosen so that its radial component A.f van- 
ishes and consequently the contribution of / A.dr along the geodesic segment from ff to 
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vanishes. It follows that BA can be expressed as 



BA\rW))- f A(r(r)).^-dr, (5) 



Eqs.©, © and (|) yield 

/dr7/{P[f(r)]exp[^{-if.fdr}] 



By inspection of Eq.(||) we arrive at: 

where Y(B) is given by Y(B) = J dffK B {Q,ff), where K B (0,ff) is the quantum ampli- 
tude for a particle of unit charge and mass to go from the origin to ff in imaginary time 
/3 in the presence of a homogenous magnetic field. This amplitude can also be expressed 



as H| 

K B (0, ff) = £ exp[-/3E n K(0K(r/), (8) 
{»} 

where u n (f) are a complete set of normalised eigenstates of the Hamiltonian H = (1/2) (— iSJ - 
A) 2 and E n are the corresponding eigenvalues. (Throughout this paper we set h = c = 1.) 
We thus arrive at the expression: Y(B) = f dff J2{n} ex P ~ \P^n]Un{0)u n (ff). 

Now we demonstrate the utility of Eq.(0) by computing the distribution of areas for 
diffusion on a plane. The function Y(B) for a particle of unit mass and charge in a 
constant magnetic field, is easily computed |l9| from the energies E n = (n + ^)B and 
the eigenfunctions u n {ff) = \f~B exp [(— B/4)rf 2 ]L n [(B/2)rf 2 ] with L n , the n th Laguerre 
polynomial, is the modulus of the vector ff. 

cosh(^-) 

From (0) we find P{B) = [cosh(/3i?/2)] _1 . Taking the Fourier transform of P{B) by 



contour integration we get the result P(A) = [/3cosh(7rj4/,9)] _1 as derived in [15]. This 
provides a check on Eq.(|7]) and illustrates its use. 

Let us now address the problem posed at the beginning of this paper: What is the 
distribution P(£l) of solid angles (fi) enclosed by a Brownian particle starting at time 
r = at the north pole r n of an unit sphere and ending at any other point ff on the 
sphere, the final point being joined to the initial one by a geodesic. (As stated earlier this 
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rule breaks down only if the final point is the south pole, which is a zero probability event). 
Unlike the planar case, P(f2) is a periodic |2(| function with period An. The generating 
function P g of the distribution of solid angles is given by 



In 



P g = / <KlP(n)e igSl (9) 



with g a half integer f2l[| . P(Q) is expressed in terms of P g by a Fourier series (rather than 
an integral) with g ranging from — oo to oo in half integer steps: 

1 oo 

P(0) = - £ e-^P g . (10) 



Relation (^) now takes the form 



47T 

9= 



where Y g is given by Y g = J df f T^ {jm} exp - [(3E 9 jtW \u* 9 m {f n )u 9 jtm {r f ) with u 9 jm {r f ) the 
normalised eigenfunctions and E 9 m the energy eigenvalues for a quantum particle of unit 
charge on a sphere subject to a magnetic field created by a monopole of quantized strength 



g [22] at the center of the sphere. The quantum numbers of the eigenstates are j, m, where 
j is the total angular momentum quantum number and m is its z component. We choose 
the vector potential in the form = g(l— cos 6), Ag = 0, which is non singular everywhere 
on the sphere except the south pole. Since Ag vanishes, (as A r did in the planar case), 
/ A along the open Brownian path does measure the solid angle as defined by the geodesic 
ruleg). 

The Hamiltonian for this problem is H 9 = -\^f e sin0& + (-^ - A 9 ^ in stan- 
dard spherical co-ordinates. The wave equation can be separated by writing ip(9, <ft) = 
exp^ipfifyRtO). From continuity of ifi, it follows that R(0) = if p^ ^ 0. Therefore, only 
those eigenstates for which p^ = 0, contribute to Y 9 . These normalised eigenstates are 
labelled by j, which ranges from \g\ to oo in integer steps. 



R 9 i(z) = f-^r^ 1 + z ) 9p ^ 9 i-^)' (i 2 ) 

where z = cos 9 and P a ' b n are the Jacobi Polynomials 23]. The energy levels of this 
system are [p^ , 23 1: E 9 j = liitk^zi .Using the integer n = j — \g\ to label the states (in 



place of j) we find 



P 9 = f)(-l)» exp(-[n(n + 1) + g(2n + l)]/9/2) 1} , (13) 
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for g > and P° = 1. From (M) we arrive at p5[ ] 



P(Q) = (1/4tt)(1 + 2 ^ <xx(Slg)P g ), (14) 

3=1/2 



where g increases in half integer steps. The function (14) is plotted numerically for various 
values of (3 in Fig. 1. The qualitative nature of these plots is easily understood. For small 
values of [3 the particle tends to make small excursions and its path encloses solid angles 
close to or Att and consequently the plots are peaked around these two values. As the 
available time (3 increases, other values of fi are also probable and the peaks tend to spread 
and the curves to flatten out. Finally in the limit of (3 — > oo the particle has enough time 
to enclose all possible solid angles with equal probability. These plots give the answer to 
the question that was raised in the beginning of the paper. 

The analytical solution given above was checked against the results of computer exper- 
iments. The diffusion process was simulated on a Silicon Graphics workstation by Monte 
Carlo methods. The simulation methodology was adapted from [|l| (where more details 
are given) and is briefly as follows. We start the simulation with 10 5 molecules located 
at the north pole of a sphere with all molecular dipoles aligned parallel to the spherical 
surface (perpendicular to the position vectors). Let Vq represent the dipoles at zero time, 
oriented along the laboratory x-axis which is perpendicular to the position vector Rq along 
the laboratory z-axis. These molecules are then subjected to random "kicks" and allowed 
to diffuse over the sphere independent of each other. The diffusion of the dipoles on the 
spherical surface is effected by rotation of the position vector R and the dipole vector V 
by the same three dimensional rotation matrix about a randomly chosen vector h normal 
to the radial position vector R. In the simulation, the diffusion constant D is chosen to be 
10 -3 per time step and the angle of rotation per time step (one iteration) is obtained from 
the corresponding probability equation for planar diffusion, which is valid for very small 
diffusion lengths. The orientation of the molecule V during diffusion is used as an index 
of the solid angle swept out by the Brownian path. We numerically evaluate this for each 
molecule applying the geodesic rule. Let Vp and Rp represent the dipole and the radial 
vectors at time (3. To determine the orientation of Vp, Vp is transformed to Vp using a 
three dimensional rotation matrix with an angle of rotation (less than n) determined by 
cos -1 (Rp .Rq) . The angle between V'p and Vq is the solid angle enclosed by the Brow- 
nian path of the diffusing molecule at time (3. Similarly, the solid angles are calculated 
for all the 10 5 dipoles at different times. These are sorted into 360 different bins, each 
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corresponding to one degree increment between and 2ir and are plotted in Figure 2. 

Since the orientation of the molecule is only determined modulo 2ir in the simulation, 
f2 is only measured modulo 2tt and not modulo 4tt. So in effect the function computed by 
the simulation is Q(f2) = P(ft) + P(0 + 2-7r). From (|l4j) it is clear that the half integer 
values of g cancel out and the integer values are doubled. Our theoretical expression for 
Q(£l) is therefore 

oo 

Q(n) = (l/27r)(l + 2^cos(n 5 )^), (15) 

9=1 

where g now increases in integer steps. The results of the numerical simulation and the 
function Q(£l) (with suitable cutoffs on g and n) given above are plotted in Figure 2 for 
different values of (3. The agreement between the simulations and the theory is excellent 
and the fractional deviation is of order 1/VN where N is the number of molecules in a 
bin. 

In computing the distribution of solid angles, we have chosen all the molecules to be 
initially at the north pole (fj = f n ) and integrated over the final point fj with a uniform 
weight, to match with the procedure followed in our simulation. In an experimental 
situation, where the diffusing particles are excited and observed with external probes, it 
may be necessary to consider a weighted average over the initial and final positions of the 
diffusing particles Q. Our analysis is easily adapted to take this into account. One just 
writes 

Y 9 = j ' dr f w f {r f )dhwi(n)K a (ri,r f ). (16) 

This leads to Y 9 = J2 n w 9 *iW 9 f exp[— f3E^/2], where w 9 n = J drw(f)u 9 l (f) and the sum is 
over all eigenstates. In this paper we have set W{(fi) = 5 2 (fi — r n ) and Wf(rf) = 1. 

In the problem of diffusion of fluorescent molecules, the tangential component of the 
orientation of the molecule determines the solid angle enclosed by the diffusing particle 
only modulo 2tt. This is because, in parallel transport of a vector, the vector rotates by 
an angle equal to the solid angle enclosed by the curve. As Pancharatnam showed Q, 
a quantum two state system transported on the Poincare sphere picks up a geometric 
phase equal to half the enclosed solid angle. This measures the solid angle modulo 4tt. 
The distribution P(£l) computed in the body of the text describes the distribution of 
geometric phases. However, since our simulations were adapted from the distribution 
measured in the simulation is Q(f2) and not P(f2). 

In [11, 14], the generating function for the distribution of solid angles for closed Brow- 
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nian paths involved only the energy eigenvalues of the Hamiltonian. In contrast, in the 
present case of open Brownian paths (closed by the geodesic rule) , one needs to use both the 
eigenvalues as well as the eigenf unctions of the Hamiltonian. This is the main difference 
between ]TT| , fbi] ] and the present theoretical analysis. 

We have computed the distribution of solid angles for Brownian motion on the sphere. 
This calculation has also been checked against computer experiments. In simple geometries 
like the surface of a sphere, an analytical treatment is possible. In more complicated 
situations, one is forced to rely on computer simulations. Monte Carlo simulations can 
be used in two distinct roles. In the present problem, where an analytical treatment is 
possible, simulations can be viewed as a computer experiment against which the theory 
can be tested. Such experiments are much easier to control than real experiments. In 
cases where an analytical treatment is not possible, the simulations serve as a substitute 
for the theory, to be checked against real experiments. The agreement between theory 
and computer experiments for Brownian motion on the sphere also serves as a check on 
the computer simulations so that these simulations can be confidently used in geometries 
which cannot be treated analytically. 
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Figure Captions 



Fig. 1: The solid angle distributions P(£l) computed from Eq. |14| at five different values 
of j3: 0.5, 1, 2, 5 and 10. The flatter curves correspond to higher j3. 

Fig. 2: Comparision between solid angle distributions Q(f2) from Monte Carlo simulations 
(red line) and the theory (Eq. |l5|) (blue line). The plot shows the distributions at five 
different values of (3: 0.5, 1, 2, 5 and 10. 
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